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Abstract 

Experiments at CEBAF will scan the intermediate-energy region of the QCD dynamics for the nucleon 
form factors and for Compton Scattering. These experiments will definitely clarify the role of resummed 
perturbation theory and of quark-hadron duahty (QCD sum rules) in this regime. With this perspective 
in mind, we review the factorization theorem of perturbative QCD for exclusive processes at intermediate 
energy scales, which embodies the transverse degrees of freedom of a parton and the Sudakov resummation 
of the corresponding large logarithms. We concentrate on the pion and proton electromagnetic form factors 
and on pion Compton scattering. New ingredients, such as the evolution of the pion wave function and the 
complete two-loop expression of the Sudakov factor, are included. The sensitivity of our predictions to the 
infrared cutoff for the Sudakov evolution is discussed. We also elaborate on QCD sum rule methods for 
Compton Scattering, which provide an alternative description of this process. We show that, by comparing 
the local duality analysis to resummed perturbation theory, it is possible to describe the transition of exclusive 
processes to perturbative QCD. 
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1 Introduction 

With the advent of CEBAF, the Continuous Electron Beam Accelerator Facility at Jefferson LAB, studies of the 
intermediate-energy QCD and analyses of various phenomena in hadronic physics, have now become a reality. 
These studies are performed experimentally with a very high beam luminosity through scattering of electrons 
on fixed targets. 

Among these, just to mention few of them, are the investigations of the hadronic resonances which are 
not classified by the "traditional" quark model, such as the J^'~' exotic states that are not accessible to any 
perturbative QCD (PQCD) prediction, and the study of the behaviour of elastic scattering at the low end of 
the intermediate-energy region. By "low-end" we refer here to the few-GeV region {Q^ « 1 — 3 GeV^), where 
PQCD does not apply, and in general, one has to resort to various models to describe the data. 

We should also mention that there is the possibility that an upgrade of CEBAF will bring up into the bulk 
of the intermediate-energy region, where the degrees of freedom of quarks and gluons become more accessible 
to a perturbative treatment, using the parton model. 
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Our aim, in this paper, is to overview, from our perspective, the phenomenology of elastic scattering off 
pions and nucleons, and to summarize the open problems, which are left for future investigations. We shall 
focus the attention on the behaviour of form factors and on Compton scattering. An experimental overview of 
these studies can be found in Ref. 

This article is entirely based on the results obtained in the last few years in ||, |[ 0, |l^, Q with the 

addition of new results, on which we elaborate in great detail. 

We remark that excellent reviews discussing the physics of form factors at intermediate energy have appeared 
recently j^, to which we refer for an overview of the subject. 

It is important to remark that the analysis of form factors at even lower energies requires the solution of 
bound state equations in the Dyson-Schwinger or Bethe-Salpeter formalism. Recent progress in the numerical 
study of form factors from the viewpoint of the Dyson-Schwinger formalism has been reported in 

In the first part of this article, we elaborate in great detail on the PQCD description of the pion and 
proton form factors, and extend the discussion to the case of Compton scattering at intermediate energy. New 
ingredients, such as the evolution of the pion wave function and the complete two-loop expression of the Sudakov 
factor, are added. We show that with the inclusion of these ingredients, the experimental data of the form factors 
at intermediate energy can be explained. 

Afterwards, we focus our discussion mainly on the connection between resummed perturbation theory and 
QCD sum rules in the case of Compton scattering. 



1.1 Comparing Two Approaches: Sum Rules versus PQCD 

The applicability of PQCD in any process requires a hard scattering contribution. In the case of the pion and 
proton wave functions, for instance, the formation of the final bound state is mediated by the hard scattering 
and we need large energy and momentum transfer for this to happen. This is the so called "asymptotic limit" , in 
which there is no direct overlap between the initial and the final state. However, there is widespread agreement 
that at current energy the description is sub-asymptotic and a pure perturbative description is not the correct 
thing to do. If this is the case, then the correct physical mechanism underlying the dynamics of elastic scattering 
in the few Gev region is the direct overlap of wave functions, which has to be described in some other way. 
One possibility to describe such "soft" behaviour is by a dispersion relation together with an Operator Product 
Expansion (quark-hadron duality). The combination of these two tools is known under the name of QCD sum 
rules. 

A preliminary local duality description of real and virtual Compton scattering was developed in Q| in the 
analysis of fixed angle Compton scattering, where a sum rule for the pion was derived. The study of the power 
corrections and of the analiticity properties of the dispersion relations for this process was presented in |^ . 
Subsequently, an analysis of Compton scattering by using resummed perturbation theory and QCD sum rules 
was proposed by two of us as a way to get information on the behaviour of the process at intermediate energy 
. It was observed that the two approaches generate different predictions for pion Compton scattering at the 
high end and at the low end of the intermediate-energy region, with an interesting interplay between them. It 
is then possible to determine the transition region to PQCD in pion Compton scattering by comparing these 
predictions. Though the analysis was limited to the pion case, the extension of this study to the proton case is 
straightforward. Work in this direction is now in progress. 

The results of indicate that sum rule predictions for the cross section of pion Compton scattering dominate 
over the resummed perturbative predictions at lower momentum transfers. At larger momentum transfers, 
instead, the perturbative predictions dominate over the sum rule ones. The analysis in [Q also showed that the 
PQCD predictions are very sensitive to variations in the momentum transfer t, while the sum rule predictions 
are less sensitive to it. 

If this interplay is indeed confirmed -as we expect- in more complex processes, such as virtual Compton 
scattering off protons (VCS), we shall have a lesson to learn: the transition to PQCD in VCS or in the real 
case (RCS) may take place at a lower scale. These conclusions, unfortunately, are affected by the usual model 
dependence from the nonperturbative parametrizations involved in both resummed perturbation theory and in 
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QCD sum rules. However, the implications are still important and can help us understand the experimental 
results. 

On the other hand, the applicability of PQCD to exclusive processes has been a matter of controversy 
for almost two decades ^. Although there is a general agreement that PQCD can successfully predict 
the behaviour of exclusive reactions as the momentum transfer goes to infinity, it remains unclear whether 
experimentally accessible energy scales (at CEBAF) are large enough to justify these predictions. In this 
work we briefly review the progress on this subject starting from the 80's, and concentrate on the resummed 
PQCD formalism developed in |ll| for elastic scattering at intermediate energy. This formalism modifies 
the conventional factorization formulas for elastic hadron scattering at large angle, which is the origin of the 
quark counting rules, by including the Sudakov effects. These effects are important both in exclusive and 
inclusive processes in special kinematic configurations in which a small transverse momentum generates large 
logarithms that need to be resummed to all orders. It was found that the effect of the resummation results 
in a suppression on the non-perturbative contributions, and improves the applicability of PQCD at the high 
end of the intermediate-energy region. It certainly remains legitimate to ask if the Sudakov effects are really 
important at, say, ~ 5 GeV^, or negligible as claimed in Q recently. 



1.2 QCD Sum Rules and Compton Scattering 

Among the semi-phenomenological methods developed for the study of the resonant region of QCD, sum rules 
are closest to a fundamental description of the strong interactions. They have been formulated in a form 
which is very similar to PQCD. In the sum rule approach a "master equation" relates the timelike region 
of external invariants in correlation functions, the region where the resonances are located, to the Euclidean 
region, where information from the Operator Product Expansion (OPE) is parameterized in terms of local 
vacuum condensates (power corrections). In practical applications only the phenomenological input from the 
lowest dimensional (quark and gluon) condensates is required. The information on the resonance region, or 
timelike region, instead, is modeled into a spectral density of states, which is described by the usual ansatz 
"lowest pole plus continuum" . 

A superconvergence condition requires the removal of the continuum contribution beyond a given "duality 
interval". In fact, equating the model spectral density to the OPE series leads to the cancellation of the 
continuum contributions, i.e., the contributions from both the left-hand side and the right-hand side of the sum 
rule The consequence of this procedure is a sum rule, which relates properties of the lowest resonance to 
an OPE series that includes both perturbative and nonperturbative terms. Radiative corrections to spectral 
functions can also be calculated 14 1. 



This way to obtain information about the resonances, which works very well for both 2- and 3-point processes, 
has largely contributed to our understanding of QCD in the intermediate-energy region In the 

extension of the method to 4-point processes, various new difficulties have been encountered, such as more 
complex dispersion relations, and the issue of the validity of the OPE. As discussed in the latter condition 
sets a constraint on the allowed range of the variables s and i, s being the center-of-mass energy, and demands 
to limit the analysis to scatterings at fixed-angle (or —t/s equivalently) . 

Once these kinematic constraints are satisfied, new information available from sum rules for 4-point processes, 
compared to 3-point ones, comes from the more general s and t dependences of the residui at double poles of 
the lowest resonance, which are related to a perturbative spectral function plus power corrections. Therefore, in 
the case of 4-point processes a true angular dependence of the sum rules comes in to play, and a parallel study 
of Compton scattering using this method [Hand the resummed PQCD factorization formalism is possible. This 
aspect has been preliminarily discussed in [kj . 

We review the above two methods that were used in the investigation of elastic scatterings, concentrating on 
the pion and proton electromagnetic form factors, and pion Compton scattering. The standard PQCD theory 
for exclusive processes is summarized in Sect. 2. The Sudakov resummation, which will be employed throughout 
the paper, is reviewed in Sect. 3, by taking the pion form factor as an example. The technique can be trivially 
extended to the proton case. The pion and proton form factors and pion Compton scattering are evaluated in 
Sects. 4, 5, and 7, respectively. We analyze pion Compton scattering using QCD sum rules in Sect. 6. Section 
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8 contains our conclusion. The Appendices contain some details of the calculations. 



2 Perturbative QCD in Exclusive Processes 
2.1 The Efremov-Radyushkin-Brodsky-Lepage Theory 

The PQCD theory for exclusive processes was first proposed by Brodsky and Lepage and by Efremov and 
Radyushkin These authors argued that an exclusive process, such as the pion electromagnetic form factor, 
can be factorized into two types of subprocesses: wave functions which carry the nonperturbative information 
of the initial and final state pions, and a hard amplitude which describes the scattering of a valence quark of the 
pion by the energetic photon. The former is not calculable in perturbation theory, and needs to be parametrized 
by a model, to be derived by nonperturbative methods, such as QCD sum rules and lattice gauge theory, or 
to be determined by experimental data. The latter, characterized by a large momentum flow, is calculable 
in perturbation theory. Combining these two types of subprocesses, predictions for the pion form factor at 
experimentally accessible energy scales can be made. 

According to this standard picture, the factorization formula for the pion form factor Ft^{Q^)^ graphically 
represented by Fig. 1(a), is written as 

Ftt{Q^)^ I dxidx2(t){x2,fi)H{xi,X2,Q,fi)(j)ixi,fi) , (!) 
Jo 

= —2Pi ■ P2 being the momentum transfer from the photon, and Pi (P2) the mometum of the incoming 
(outgoing) pion. /i is a renormalization and factorization scale, below which the QCD dynamics is absorbed 
into the pion wave function (j), and above which the dynamics is regarded as perturbative and absorbed into the 
hard amplitude H . The pion wave function <^(a;, n) gives the probability of a valence quark carrying a fractional 
momentum xP in the parton model at the energy scale fj,. Usually, the ^ dependence of (j) is neglected. H 
is obtained by computing the quark-photon scattering diagrams as shown in Fig. 1(b). All the lines in (j) are 
thought of as near or on the mass shell, while those in H are far off-shell by the scale Q. 

To make predictions for F^r, one substitutes the lowest-order expression of H obtained from Fig. 1(b), 

H{xi,X2,Q,l-i.) = 7^5—' 2 

XiX2Q^ 

asipi) being the running coupling constant, and the asymptotic pion wave function [ pO| , 

0^^(x) = ^a;(l-x), (3) 



into Eq. (|), where Cf = 4/3 is a color factor, Nc — i the number of colors, and /^r — 0.133 GeV the pion 
decay constant. Because of the behaviour of (j)^^ , the main contributions to Eq. (Q) come from the regions 
with intermediate xi and X2. Higher-order corrections to H then produce the logarithms of the type ln((5^//i^), 
which may be so large as to spoil the perturbative expansion. To eliminate these logarithms, a natural choice 
of is ^ q2 jg gg^gy obtain Q'^F^ - 0.12 GeV^ for _ 4_io GeV^, which is only 1/3 of the data 
0.35 GeV^ ^ This contradiction implies that the pion form factor has not yet become completely 

asymptotic at experimentally accessible energy scales. 

Hence, one may resort to a preasymptotic model, the Chernyak and Zhitnitsky (CZ) wave function derived 
from QCD sum rules pst , 

cj^^^{x)^^x{l-x){l^2xf , (4) 

which possesses maxima at the end points a; — > and x — » 1. It was found that the CZ wave function enhances 
the predictions for Q^P^ to 0.30 GeV^, and improves the match with the data. 
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In spite of the success of the standard PQCD theory, the perturbative expressions for the pion form factor 
were brought under searching criticism by Isgur and Llewellyn Smith and by Radyushkin [ p6| . Since cj)'-^^ 
emphasizes the contributions from small xi and X2, the end-point regions are important. They argued that 
higher-order corrections to H in fact produce logarithms like \ii{xiX2Q^ / ^J?) and thus the natural choice of 
fj? is fi^ — XiX2Q^- For intermediate Xi and X2, this choice is consistent with ji^ = employed above. 
For small Xi and X2, it results in as{xiX2Q'^) > 1, indicating that the end-point regions are nonperturbative 
regions. A careful analysis reveals that the enhancement of the predictions for Q^iv by 0*-^^ is due to the 
amplification of the nonperturbatve end-point contributions. As a consequence, the perturbative calculation 
loses its self-consistency as a weak-coupling expansion. 

A cutoff to avoid the divergences of was proposed : 

TT 

" /?oln[(A*2+4m2)/A2] ' 

with A = Aqcd the QCD scale parameter and (3q = (33 — 2n/)/12, n/ = 3 being the number of quark 
flavors. rUg is interpreted as a dynamical gluon mass, which is acquired from the low momentum region of 
those radiative corrections that induce the running of The leading-order PQCD predictions are then 

stabilized at low momentum transfers because of the frozen In this interesting approach the value of nig 
has to be determined by matching the predictions with the data. However, the results from Eq. (0) are slightly 
sensitive to the variation of nig around the best choice « 0.3 GeV^ [|2^, and thus the theory with this cutoff 
becomes less predictive. An alternative way to freeze the coupling constant at low momentum transfer is the 
BLM method, in which the renormalization scale 11 is fixed according to some criteria. We refer to for more 
details on the application of this method to exclusive processes. 



2.2 Resummed PQCD Formalism 

The above discussions hint that the end-point regions should be treated in a different way. In general, a valence 
quark in the pion can carry a small amount of transverse momenta kx- For intermediate xi and X2, fcr that flow 
from the wave functions through the hard scattering are neglected as power-suppressed corrections, as shown 
by the hard gluon propagator 

XiX2Q'^ + k'^ XiX2Q^ \ XiX2Q^ 

The first term is exactly the asymptotic hard scattering in Eq. (^, and the second is suppressed by powers 
of at fixed x. If the end-point regions are not important, one may drop the second term. However, the 
end-point difficulties indicate that the higher-power effects are crucial, and should be kept at the outset in the 
derivation of the new perturbative expression for the pion form factor. 

The introduction of the kx dependence leads to three observations immediately. First, the pion form factor 
becomes a two-scale (Q and kx) problem, and the resummation technique may be required to organize the 
large logarithms hi{Q/kT) from radiative corrections contained in the wave function. Second, we must analyze 
the process in the Fourier transform space of kr, denoted by b, which is regarded as the transverse separation 
between the valence quarks of the pion. Third, the quantity 1/6, now considered as one of the characteristic 
scales of the hard scattering amplitude, should be substituted for the argument of as if l/b^ > xiX2Q^ ■ The 
nonperturbative regions are then characterized by the regions of large b and small x. 

We find that the resummation of the large logarithms behaves as 

exp[— const, x In Q ln(ln Q/ ln6)] , (7) 

which suppresses the elastic scattering at large spatial separation. This property, called Sudakov suppression 
[ p^ |29| , |30| , makes the nonperturbative contributions from large 5, no matter what x is, less important, without 
introducing any phenomenological parameters such as a gluon mass m,g. The relevance of the Sudakov effects 
was pointed out by Lepage and Brodsky in . The underlying physical principle is that the elastic scattering of 
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an isolated colored parton, such as a quark, is suppressed at high energy by radiative corrections. We shall show 
that the new perturbative expression reduces to the standard one in Eq. (|l]) as Q ~* oo? but, at lower momentum 
transfers, it takes into account the infinite summation of those higher-order effects whose b dependence makes 
the perturbative theory more self-consistent -at least at the high end of the experimentally accessible energies-. 
This self-consistency means that numerical results come predominantly from the momentum regions in which 
the effective coupling is relatively small. The above argument can be quantified by studying the dependence of 
the contribution to the form factor on a cutoff in 5, instead of in x |||, ^ , from which the importance of the 
perturbative region is examined. It is observed that the inclusion of Sudakov corrections produces numerical 
effects very similar to the approaches with a frozen as at low energies, and in a good agreement with the data. 
A brief review of the standard asymptotic expression for the form factors is referred to p7|, |3^ , |3^ . 

The perturbative calculations should be thought of as complementing the QCD sum rule approach [H , which 
is one of the topics of this paper. 



3 Sudakov Resummation 

As we have discussed above, the applicability of PQCD to exclusive processes at intermediate energy can be 
improved by including the transverse momentum dependence and the Sudakov suppression into the factorization 
formulas. In this section we review the resummation technique, taking the pion form factor as an example. We 
investigate radiative corrections to the basic diagrams for the form factor, and explain how they are factorized. 
For details on the Sudakov resummation, we refer to pg] . 



3.1 Factorization in b Space 

To construct the factorization of QCD processes we first isolate the leading momentum regions of those radiative 
corrections from which important contributions to the loop integrals arise. There are two types of important 
contributions: coUinear, when the loop momentum is parallel to the pion momentum, and soft, when the loop 
momentum is much smaller than the momentum transfer Q^. A small amount of transverse momenta kx is 
associated with the valence partons that enter hard scattering which is taken as an infrared cutoff. Each type of 
the important (soft or coUinear) contributions gives large single logarithms, and their overlap generates double 
(leading) logarithms, which are a characteristic feature of gauge theories. These large logarithms, appearing in 
a product with as, must be organized in order not to spoil the perturbative expansion. Single logarithms can 
be summed to all orders using rcnormalization group (RG) methods, while double logarithms must be treated 
by the resummation technique IH] . 

The diagrams shown in Fig. 2 represent 0(q!s) corrections to the basic factorization of the pion form factor, 
containing the large logarithms mentioned above. In the axial (physical) gauge n ■ A = 0, n being a gauge 
vector and A the gauge field, the two-particle reducible diagrams in the channels of the external pions, like 
Figs. 2(a) and 2(b), have double logarithms from the overlap of coUinear and soft divergences, while the two- 
particle irreducible corrections, like Figs. 2(c) and 2(d), contain only single soft logarithms. This distinction 
is consistent with the physical picture of a hard scattering process: two partons moving in the same direction 
can interact with each other through coUinear or soft gluons, while those moving apart from each other can 
interact only through soft gluons. The Sudakov corrections come entirely from diagrams that are two-particle 
reducible. Below we shall concentrate on reducible corrections, and demonstrate how they are summed into a 
Sudakov factor. 

A careful analysis shows that soft divergences cancel between Figs. 2(a) and 2(b), as well as between 2(c) 
and 2(d), in the asymptotic region with the transverse separation 6 — > 0. Therefore, reducible corrections 
are dominated by coUinear divergences, and can be absorbed into the pion wave funtion, which involves a 
similar dynamics. Irreducible corrections, due to the cancellation of soft divergences, are then absorbed into the 
hard scattering amplitude H. Hence, the factorization picture holds at least asymptotically after the radiative 
corrections are included. The cancellation of soft divergences is closely related to the universality of wave 
functions. For a large 6, double logarithms are present and the resummation technique must be implemented. 
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When the transverse degrees of freedom of a parton are included, the factorization of the reducible corrections 
into the pion wave function must be performed in b space. Fig. 2(a), where the loop momentum I does not flow 
through the pion wave function, is computed directly. For Fig. 2(b), where I may be routed through the wave 
function tp, the associated integrand contains a factor ^p{xi + l^/P^ , \kiT + 1t|)- Here we have assumed that 
l~ is routed through the outgoing pion wave function. Since the behavior of ip is unknown, the loop integral 
can not be done directly. We employ the approximation ip{xi + l~^/P^ , \hiT + 1t|) ~ tpi^i, |kir + l^j), since 
P^ is large and the dependence is negligible. However, the transverse momentum kiT is not large, and thus 
the It dependence remains. This difficulty can be overcome by introducing a Fourier transform, 

^^expH(kiT + lT)-b]7'(a;i,6) , (8) 

with V the Fourier transformed wave funciton. The factor exp(— ily • b) is absorbed into the loop integrand, 
which then becomes computable, and the factor exp(— zkiT^ • b) Fourier transforms the hard part H to b space. 
Using the above reasoning, the factorization formula for the pion form factor is written as 

/•I r (^2^5 

F-,t{Q^) = J dxidx2 J j^^V{x2,b,P2,fi) 

xH{xi,X2,b,Q,fi)'P{xi,b,Pi,fi), (9) 

where the expression of the Fourier transformed hard amplitude H will be given in next section. Equation (^) 
depends only on a single parameter 6, because the virtual quark line involved in H is thought of as being far 
from mass shell, and shrunk to a point The wave function V is defined in terms of explicit matrix elements, 
which complement the diagrammatic descriptions presented in Fig. 1(a) [ p^ , 

Vix,b,P,^i) = J d^kTe'^^-^i:{x,kT,P,fi) 

^e"^+^"(0|r(g-(0)7+9(y-,0+,b))|^(P)) . (10) 

V includes all leading logarithmic enhancements at large 6, which will be resummed below. 

Equation (||) is an intermediate step in deriving Eq. (|l|) |Q. If one sets b to zero in V, and integrates over 
b in H, the former reduces to the latter. For fixed Xi ^ 0, the two expressions are equivalent up to corrections 
that fall off as a power of Q. However, Eq. (^) retains much more higher-power information at the limit — > 
for fixed Q, which is potentially important. 

3.2 Technique 

The basic idea of the resummation technique is as follows. If the double logarithms appear in an exponential 
form in Eq. (^), the task will be simplified by studying the derivative dV/d\nQ = CV. The coefficient function 
C contains only large single logarithms, and can be treated by RG methods. Therefore, working with C one 
simplifies the double-logarithm problem into a single-logarithm problem. 

The two invariants appearing in V are P ■ n and n^. Due to the structure of the gluon propagator in the 
axial gauge. 



n - I (n-iy 

V depends only on a single large scale [P ■ n)^ jv? . We choose the Breit frame such that P^ = P2 = Q/V^ 
and all other components of P's vanish. It is then easy to show that the differential operator d/d\nQ can be 
replaced by d/dn, 

d d 

-T' = --^P^-^T' . (12) 



dlnQ P-n dn° 
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The motivation for this replacement is that the momentum P flows through both quark and gluon hnes, but n 
appears only in gluon lines. The analysis then becomes simpler by studying the n, instead of P, dependence. 
Applying d/dua to the gluon propagator, we obtain 

-iLAT/^'' = _J_(7V^"r + 7^^°?^) . (13) 
dua I ■ n 

The momentum (/'') is contracted with the vertex the "differentiated" gluon attaches, which is then replaced 
by a new vertex, 

gT'^-^f—P^ , (14) 
P ■ ni ■ n 

T° being a color matrix. This new vertex can be easily read off from the combination of Eqs. ( |l2| ) and (p^). 
After adding together all the diagrams with "differentiated" gluons and using the Ward identity, we arrive at the 
differential equation of V shown in Fig. 3(a), in which the new vertex is represented by a square. An important 
feature of the new vertex is that the gluon momentum I does not lead to coUinear divergences because of the 
nonvanishing n^. The leading regions of I are then soft and ultraviolet, in which Fig. 3(a) can be factorized 
according to Fig. 3(b) at lowest order. The part on the left-hand side of the dashed line is exactly V, and that 
on the right-hand side is assigned to the coefficient function C. 

Hence, we introduce a function /C to absorb the soft divergences from the first two diagrams in Fig. 3(b), and 
a function G to absorb the ultraviolet divergences from the other two diagrams. The soft subtraction employed 
in Q ensures that the involved momentum fiow is hard. Generalizing the two functions to all orders, we derive 
the differential equation of V, 

^^7'(x,6,P,Ai) - [2IC{bfi)+g{xiy/fi)+g{{l^x),y/fi)] 

xV{x,b,P,^i) . (15) 

/C and Q have been calculated to one loop, and their single logarithms have been organized to give the evolutions 
in b and Q, respectively [ p^ . They have ultraviolet poles individually, but their sum IC + G is finite. 
Substituting the expressions for IC and G into Eq. (^5|), we obtain the solution 

V{x, b, P, fi) = exp [-s{xQ, 1/b) - s((l - x)Q, 1/b)] P{x, b, y) . (16) 

The exponent s{(,Q, 1/b) for S, — x and 1 — a; is written as pO| | 

A{as{p)) + B{as{p)) 



s{^Q,l/b) 



Ji/b P I \V2pJ 



(17) 



n/b 

where the anomalous dimensions A to two loops and B to one loop are given by 



A = 



Cf— + 

TT 


'67 


^2 


10 


."9" ^ 


3 


Uf 

27 ^ 


^^Inf 

3 TT V 




-) 




2 




1 



2 

B = ^-ln( ) , (18) 



with 7b the Euler constant. The two-loop expression of the running coupling constant, 

aa(M) ^ 1 Pi lnln(A^VA2) 

IT ~ /?oln(/.VA2) Pi \n\^P/K^) ' 

will be substituted into Eq. (0), with the coefficients 



(19) 



A = = ,20) 



To derive Eq. (|17D, a space- like gauge vector n oc (1,-1,0) has been chosen. We require the ordering of the 
relevant scales to be > 1/6 > A as indicated by the integration of p from 1/6 to in Eq. (p^). The QCD 
dynamics below 1/6 is regarded as being nonperturbative, and absorbed into the initial condition 'P(^, 6, ^). 
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4 The Perturbative Pion Form Factor 

In this section we evaluate the pion form factor in the framework developed above. Compared to the analysis 
in | [Tl| , the more complete Sudakov factor derived from the two- loop running coupling constant as in Eq. ( p^ ) 
will be employed and the evolution of the pion wave function will be taken into account. The above resummed 
PQCD formalism has been widely applied to various exclusive processes, such as the time-like pion form factor 
|p4[ , proton-antiproton annihilation [pSl and proton-proton Landshoff scatterings pq| . 



4.1 Evolution in b 



We continue the organization of the large logarithms contained in the factorization formula for the pion form 
factor in Eq. (||). The functions V and H still contain single logarithms coming from ultraviolet divergences, 



which need to be summed using the RG equations |10 
VH{xi, b, Q, n) = A-iqH{xi,b, Q, ^) , 

with 



(21) 
(22) 



(23) 



and 7^ = — a^/Tr being the quark anomalous dimension in the axial gauge. Solving Eq. (pl|), the large-6 behavior 
of V is summarized as 



V{x, b, P, Ai) 



exp 



-s{xQ,\/b) 



s{{l-x)QMb)^2 r ^7,(«.(m)) 
Ji/b M 



xVix,b,l/b) , 



(24) 



where the arguments b and l/b of the initial condition V denote the intrinsic and perturbative b evolutions, 
respectively. Below we shall ignore the intrinsic b dependence, and assume that V = (j){x, l/b). 
The RG solution of Eq. 0) is given by 



H{xi,b,Q, fi) = exp 



-4 



H{xi,b,Q,t) , 



where t is the largest mass scale involved in the hard scattering. 



t — max (-^x 1X2(3, l/b) 



(25) 



(26) 



The scale yJx\XiQ is associated with the longitudinal momentum of the hard gluon and l/b with the transverse 
momentum. The anomalous dimensions 7^ allow us to take into account the strong couplings of quarks when t 
is small. With the large logarithms organized, the initial condition H{xi, b, Q, t) can be taken to coincide with 
its lowest-order expression, obtained by a Fourier transform of the hard scattering amplitude H in momentum 
space, as derived from Fig. 1(b), 



H = 



WnasCpxiQ^ 



(xiQ^ -I- ]s.Ij.){xiX2Q^ 

IGnasCp 

XiX2Q^ + (kiT - k2T)' 



k2T)2) 



(27) 
(28) 



We have neglected the transverse momentum in the numerator of the first expression, and neglected, in addition, 
the transverse momentum carried by the virtual quarks in the denominator of the second expression. These 
quark propagators are linear rather than quadratic in x. Because Eq. (psj) depends on the combination of the 
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transverse momenta, the factorization formula for the pion form factor involves only a singlc-& integral as shown 
in Eq. (||). Note that Eq. ( p8|) coincides with Eq. (^) at k^T = 0. 
Inserting Eqs. and (|25|) into Eq. (||), it becomes 



IBttCf / dxidx2(t){xi,l/b)(l){x2,l/b) / bdbas (t)Ko{^Xi X2Qb) 
Jo Jo 

X e-Kp[-S{xi,X2,b,Q)] , 



(29) 



with the complete Sudakov exponent 



S{xi,X2,b, Q) = V [s(x,Q, 1/b) + s((l - xOQ, l/b)] + 4 / ^7,(0,(^2)) 



(30) 



i^o is the modified Bessel function of order zero, which is the Fourier transform to b space of the gluon propagator. 
Equation ( p^ ) for the pion form factor, being a physical quantity, is explicitly /^-independent. This is the 
consequence of the RG analysis in Eqs. ( pl| ) and (p^), and the advantage of the resummed PQCD formalism. 

Another advantage of Eq. (|2|), compared to Eq. (Q), is the extra b dependence. For x close to zero, i.e., 
xiX2Q^ < 1/^^, the scale of the radiative corrections is dominated by the transverse distance bridged by the 
exchanged gluon. If this distance is small, radiative corrections with the argument of set to t in the hard 
scattering amplitude will be small, regardless of the values of x. Of course, when b is large and xiX2Q^ is small, 
radiative corrections are still large. However, the Sudakov factor exp(— 5), exhibiting a strong fall off at large 
b, vanishes for b > 1/A. We expect that, since Eq. (Q) is the asymptotic limit of Eq. (|29|), as Q increases, 
the b integral will become more and more dominated by the small-6 contributions. If the Sudakov suppression 
is strong enough, such that the main contributions to the factorization formula come from the short-distance 
region, the perturbative expansion will be rendered relatively self-consistent. Hence, we propose to analyze 
Eq. (p9|), not by cutting off the x integrals near their endpoints, but by testing its sensitivity to the large-& 
region. 

The evolution of the wave function </> with b is written as [p8|. 



0(x,l/6) = -^i^x{l-x) 



Y,ar,Cl'\2x-l) 



.n=0 



as{l/b) 

Olsip-o) 



with /xo ~ 0.5 GeV. The coefficient and the exponent 7„ are 
5(2n-f 3) 



(n-t-2)(n-|- 1) 



1 



33 - 2n 



f 



n+l 

^ k (n+l 



k=2 



(n+ l)(n + 2) 



with /o = 2/15, Ii = 0, I2 = 8/35 and /„ = for n > 3. The relevant Gegenbauer polynomials are 



C',/'ix)^l, C'^'^ix)^-yx ^ 
That is, we adopt the pion wave function 



-.3/2. 
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(/.(X, l/b) 



3/. 



:X{1 — X) 



(31) 

(32) 
(33) 

(34) 

(35) 



It is easy to observe that Eq. ( ^5| ) approaches <^^^ in Eq. (||) as 1/6 ^ 00, and 0*^^ in Eq. (^) as 1/6 —> /ip. 
This is the reason (\)^^ is an asympotic wave function, while t/)'^'^ is a preasympotic one. 
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We are also interested in the full expression of H in Eq. (27), except for the approximate one in Eq. (|28| ) 



The corresponding factorization formula involves a double-6 integral 

pl poo 

Fti{Q^) = IGttCf / dxidx2 / bidbib2db2<p{xi,l/bi)(f>{x2,l/b2) 
Jo Jo 

xas{t)Ko{y/xiX2Qbi) exp 0:2, 61, 62, Q)] 

X [9{b, - b2)Ko{V^Qb,)Io{V^Qb2) + 

e{b2-b,)Ko{^Qb2)Io{V^Qbi)] , (36) 

with the Sudakov exponent 
2 

S = ^[s(a:,Q, 1/6,) + s((l-x,)Q, 1/60] 

i=l 

+2/* ^j^{asifl)) + 2 f ^j^ia.ifi)), (37) 
and the hard scale 

i = max [Vil^Q, 1/61, 1/62] . (38) 
4.2 Numerical Results 



Results on the behaviour of the Sudakov suppression in the large 6 region where presented in |11|. There is 
no suppression for small 6, where the two fermion lines are close to each other. In this region the higher- 
order corrections should be absorbed into the hard amplitude, instead of being adsorbed into the pion wave 
function. Hence, we set any factor exp[— s(^(5, 1/6)] to unity for ^(3/\/2 < 1/6. Similarly, we also set the 
complete Sudakov exponential exp(— 5) in Eq. ( |29| ) to unity, whenever it exceeds unity in the small-6 region. 
As 6 increases, exp(— S*) decreases, reaching zero at 6 = 1/A. Suppression in the large 6 region is stronger for 
larger Q. 

We choose the QCD scale A = 0.15 GeV for the numerical analysis of the pion form factor below. It can be 
shown that the results will differ only by 10%, if A = 0.1 GeV or A = 0.2 GeV is adopted. 

In order to see how the contributions to Eq. ( [29| ) are distributed in the 6 space under Sudakov suppression, 
the integration is done with 6 cut off at be- Typical numerical results arc displayed in Fig. 4 for the use of in 
Eq. (|35|). The curves, showing the dependence of Q^Fj^ on 6c, rise from zero at be = 0, and reach their full height 
at be = 1/A, beyond which any remaining contributions are regarded as being truly nonperturbative. To be 
quantitive, we determine a cutoff up to which half of the whole contribution has been accumulated. For — 4 
and 9 GeV^, 50% of Q'^F^ comes from the regions with 6 < 3.6 GeV^^ [as{l/b)/Tr < 0.32] and with 6 < 3.0 
GeV~^ [as(l/6)/7r < 0.21], respectively. A liberal standard to judge the self-consistency of the perturbative 
method would be that half of the results arise from the region where as{l/b)/TT is not larger than, say, 0.5. 
From this point of view, PQCD becomes self-consistent at about ~ 4 GeV^. This conclusion is similar to 
that drawn in |ll| , where (f>^'^ and (j)"^^ were employed. As anticipated, the perturbative calculation improves 
with an increase of the momentum transfer. 

Results of Q^Ft^ from c/)^^ , (j>'~'^ and cj) in Eq. f^), as well as the experimental data |2^, are 
shown in Fig. 5. The curve corresponding to Eq. ( |35[ ) is located beyond the one corresponding to (j)^^ (and 
beyond the one corresponding to (j)^^ , of course), indicating that the evolution of the pion wave function, 
especially the portion from 1/6 < /iq, plays an essential role. Our predictions for > 2 GeV^ are in a good 
agreement with the data. For < 2 GeV^, the curve deviates from the data, implying the dominance of 
the nonperturbative contrubution. However, we emphasize that the wave function (j)(x,l/b), regarded as an 
expansion in as{l/b)/as{p-o), becomes meaningless gradually as 1/6 ^ A < /ip. Hence, we shall not claim that 
the match of our results with the data is conclusive. It is possible that higher-twist contributions from other 
Fock states are important as argued in |3^. Nonetheless, the resummed PQCD predictions are indeed reliable 
at least near the high end of the accessible energy range, where the curve from <f> approaches that from cp'-'^ . 
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At last, we evaluate the pion form factor using the double-6 factorization formula in Eq. (36), and the results 
are presented in Fig. 6. Those derived from the single-6 formula in Eq. (^9|) and the experimental data are also 
displayed for comparison. It is found that the two curves are consistent with each other especially in the region 
with ^ 5-10 GeV^. The consistency hints that the approximation of neglecting the transverse momenta 
carried by the virtual quarks makes sense. Therefore, we shall employ the same approximation in order to 
simplify the multidimensional integrals, when evaluating the hard scattering amplitudes for the proton form 
factor. 

We have stated that the function Vix, b, l/b) in Eq. ( |24| ) contains the intrinsic b dependence denoted by its 
argument b, except for the perturbative evolution with b shown in Eq. (^5|). This intrinsic dependence, providing 
further suppression at large b, has been discussed in detail in ]40| . Since it depends on the parametrization, we 
shall not consider it here. 

5 The Perturbative Proton Form Factor 

With the success of the resummed PQCD formalism in the study of the pion form factor, we are encouraged 
to apply it to the proton form factor. We concentrate on the proton Dirac form factor Ff , which has been 
widely calculated based on the standard PQCD theory §, |8[ ||, ||^, || |l|, ||, H, H- Another form factor 
-p2(Q^) ~ 0{1/Q^) is much smaller than Fi{Q^) ~ 0(1/(3'*) in the asymptotic region |25 . It has been shown 
that the asymptotic proton wave function (j){xi) cx X1X2X3 fails to give the correct sign for Ff, and generates 
much smaller values than the experimental data Q. Highly asymmetric distribution amplitudes, such as 
the Chernyak-Zhitnitsky (CZ) and King-Sachrajda (KS) |Q wave functions derived from QCD sum rules, 
reverse the sign and enlarge the magnitude of the predictions. However, it was also critically pointed out that 
the perturbative calculations based on the above phenomenologically acceptable wave functions are typically 
dominated by soft virtual gluon exchanges, and violate the assumption that the momentum transfer proceeds 
perturbatively We shall show that the resummed PQCD formalism can give reliable predictions which are 
consistent with the experimental data. 

5.1 Factorization Formula 

As stated in Sect. 3, the factorization formulas for exclusive QCD processes require the investigation of those 
leading regions in the diagrams describing the radiative corrections where the dominant contributions arise 
from. Similarly, we work in the axial gauge n • ^ = as in the pion case. The two-particle reducible diagrams, 
like Figs. 7(a) and 7(b), which are dominated by the coUinear divergences, are grouped into the distribution 
amplitudes, and generate their Sudakov evolution. The bubbles may be regarded as the basic photon-quark 
scattering diagrams. The two-particle irreducible diagrams Figs. 7(c) and 7(d) contain only soft divergences, 
which cancel asymptotically. These types of corrections are then dominated by momenta of order Q^, and 
absorbed into the hard scattering amplitude H. Hence, the factorization formula for the proton Dirac form 
factor includes three factors Q 

FfiQ^) = I (dx)(dT')(dkT)(dk^)Fo'/3'y(fc,',P',Ai) 
Jo 

xHo,- p''^-aP'i{ki,K^Q, y)Yap^{ki, P, ^j.) , (39) 

with P — (P+,0,0) the initial-state proton momentum, and 

3 

[dx) = dxidx2dx'i5(^y^^Xi — 1) , 

i=l 

3 

(dkr) = d2kiTd^k2Td^k3T(5(X] ^^t) • (40) 
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Xi — kf / and k^T are the longitudinal momentum fractions and transverse momenta of the parton j, 
respectively. The primed variables P' — (0, P^,0) and x[ — k~/P~ are associated with the final-state 
proton. = 2P ■ P' is the momumtum transfer, and the scale ^ is the renormalization and factorization scale. 
In the Breit frame we have P+ =P'- = Q/V2. Similarly, we have taken into account the transverse decrees of 
freedom of a parton as indicated in Eq. (p9[). If the end-point regions are not important, one may drop the kx 
dependence in the hard scattering H , and integrate out kx in the proton distribution amplitude Y, arriving at 
the standard factorization formula. 

The initial distribution amplitude Yap^, defined by the matrix element of three local operators in the axial 



gauge, is written as ^ 



2 

[{VC)Ml5N)^V{h,P,fi) + {lf^5CUN^A(h,P,fi) 



- {<y^..P''C)^pij^j5N)^Tih, P, fi)] , (41) 

where |P) is the initial proton state, u and d the quark fields, a, b and c the color indices, and a, [3 and 7 the 
spinor indices. In our notation, 1,2 label the two u quarks and 3 labels the d quark. The second form, showing 
the explicit Dirac matrix structure, is obtained from the spin property of quark fields psf , where N is the proton 
spinor, C the charge conjugation matrix and cr^^ = [^^,^^\/2. The dimensional constant /at is determined by 
the normahzation of the distribution amplitude [^. The amplitude Fa'/3'7'(fc^, P', ^) for the final-state proton 
is defined similarly. The permutation symmetry between the two u quarks gives the relations among the wave 
functions 

V{ki,k2,k3,P,^) = F(fc2,fci,/c3,P,^) , 

A{ki,k2,k3,P,fj,) = -A{k2,ki,k3,P,^j,) , 

T{ki,k2,k3,P,fi)=T{k2,ki,k3,P,fi) . (42) 
Since the total isospin of the three quarks is equal to 1/2, we have 

2T{ki,k2,k3,P,^j,) = ip{ki,k3,k2,P,fJ,) + ■ip{k2,k3,ki, P, fj.) , (43) 

with 

iP{ki,k2,k3,P,fi) V{ki,k2,k3,P,n) - A{ki,k2,k3,P,fi) . (44) 
The combination of Eqs. (M2) and (B3h leads to 



V{ki,k2,k3,P,fi) = i [^{k2,ki,k3,P,fi) + -ipiki, k2, k^, P, fi)] 

A{ki,k2,k3,P,n) = ^ ['4!{k2,ki,k3,P,fi)-tl;{ki,k2,k3,P,n)] . (45) 

Hence, the factorization formula in fact involves only one independent wave function ip. 

The hard-scattering kernel Ha'p'-y'ap-y can be calculated perturbatively. To lowest order in with two 
hard exchanged gluons, 42 diagrams are drawn for the proton form factor. This number is reduced to 21 when 
the permutation symmetry between the incoming and outgoing protons is considered. These 21 diagrams are 
further divided into two categories, which transform into one anoother by an interchange of the two u quarks, 
so that it is enough to calculate only 11 of them, as shown in Fig. 8. The explicit expression of the integrand 
Ya' aw Ha' fS'^'a/s-yYafs-y, Corresponding to each diagram in Fig. 8, is listed in Table 1. For more details we refer 
to [He . The dia grams with three-gluon vertices vanish at leading twist and are thus ignored. 
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Applying a series of variable changes |46|| , the summation of the contributions over the 42 diagrams is reduced 
to only two terms. Equation (^) becomes 

F?iQ') = ^ l\dx){dx'){dkr){dk'^)[fN{f^)r 

2 

X ^ i/,(fc„fc^,Q,Ai)*,(fc„fc:,P,P',/i) . (46) 



with 



Hi 



3[(1 - Xi){l - x'i)Q2 + (kiT - k;^)2][a;24Q2 + (J^aT " KtY] ' 



= '^-^ (47) 



and 



5- 



2(^V'')l23 + 8(TT')l23 + 2(V'V'')l32 + 8(rr')l32 ^ (V'V'')321 " (V'V'O 



231 



(l-xi)(l-a;;) 



^ 2(V;V^0l32 - 2(rr)i23 (V^Vy)l23 - 8(rT0i32 - 2(7/^^^321 . . 

{l-x^){l-x',) ^ (i_a;3)(l-a;'i) • ^ 

The notation {ipijj') 123 represents 

(V'V'')i23 = ^(fci, fe, A:3, P, fc2, ^3: : (49) 

and (TT') 123 is defined in a similar way and we have inserted the electric charges of the quarks. The transverse 
momenta carried by the virtual quarks in the hard scattering subdiagram have been neglected for simplicity as 
indicated in Eq. (^). This approximation equates the initial and final transverse distances among all pairs of 
valence quarks. At the same time, fermion energies such as XiQ^ will not be among the characteristic scales of 
the hard scattering, and the evolution of iJ in Q will be simpler. The accuracy of this approximation has been 
roughly discussed in Sect. 4. 



5.2 Sudakov Suppression 

Reexpressing Eq. ( p6| ) in the Fourier transform space, we obtain 

^f(Q') = ^(dx){dx'){dh)[f^{^,)f 

X Hj {xi , X ■ , b,,, Q, {xi , X ■ , bj , P, P', ^) , 



(50) 



with hi the conjugate variables to k^T, and (dh) — c?^bi(i^b2/(27r)^. The above factorization formula involves 
only the integration over b (not b') as stated above. The Sudakov resummation of the large logarithms in the 
Fourier transformed wave function V leads to 



'P{xi,hi,P,fi) 



exp 



s{xiQ,w) - 3 



1=1 



{as{^J■)) 



X(j){Xi,w) , 



with the scale 

w = min(l/6i, I/62, l/^^s) , 



(51) 
(52) 
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and 63 = |bi — b2|. The function 0, obtained by factoring the Q dependence out of corresponds to the naive 
parton model. Again, we have dropped the intrinsic 6 dependence of the proton wave function. The exponent 
s{xQ,w) has been given in Eq. (p7|). 

Note the choice of the infrared cutoff in E q. (p^ ) , which differs from that adopted in the previous analysis 
but the same as w — 1/ max(6i, 62, 63) in [[f8|. In |46| the different transverse separations were assigned to 
each exponent s and to each integral involving 7^: 



'P(xi,bi,P,/i) 



exp 

L ;=i 

X(/)(Xj, w) 



^ Ls(x,Q,l/5;) 



i/b, M , 



(53) 



The Sudakov factor in Eq. (|5^) does not suppress the soft divergences coming from 6; 1/A completely. For 
example, the divergences from 61 — > 1/A, which appear in the integral involving 7^ and in (j){xi,w) at w — > A, 
survive as xi —^ 0, since s{xiQ, I/61) vanishes, while s{x2Q, 1/^2) and s^x^Q, I/63) may remain finite. On the 
other hand, w should play the role of the factorization scale, above which QCD corrections give the perturbative 
evolution of V shown in Eq. (|5^), and below which they are absorbed into the nonperturbative initial condition 
(j). Hence, it is not appropriate to choose the infrared cutoff of the Sudakov evolution different from w. An 
alternative choice of the scale ISStl, 



w' = max(l/6i, I/62, 1/63) , 



(54) 



also removes the soft divergences. We argue that the factorization scale should be chosen to diminish the largest 
logarithms contained in (f), which are proportional to ln(/i/w). However, w' does not serve this purpose. In 
the numerical analysis below we shall truncate the integrations over bi and 62, in order to examine how much 
contributions are accumulated in the short-distance (perturbative) region. Since small values of bi and 62 do not 
guarantee a small 63 due to the vector sum, we shall require 1/63 to be greater than A, such that the running 
coupling constant as{w) is well defined. 

The evolution of the hard scattering amplitudes is written as 



Hj{xt,Xi,hi,Q, fi) = exp 



1=1 ^ 
xHj{xi, x\^ b,j, Q, tji, tjq) 



(55) 



with 



^{l~xi){l~x'^)QMbi 



\/xix\Q, l/bi 



ii2 = ^22 = niax 



\fx^^Q, l/b-2 



(56) 



The first scales in the brackets are associated with the longitudinal momenta of the hard gluons and the second 
with the transverse momenta. The arguments tji and tj2 of Hj mean that each as is evaluated at the largest 
mass scale of the corresponding hard gluon. 

Inserting Eqs. (^ and ( |55| ) into (|50|), we obtain 



2 A /•! I-OD 

,=1 ^' "'0 Jo 



2ir 

db2 I de[fN{w)f 




xHj{xi,x'i,bi, Q, tji,tj2) '^]{xi, x\, w) 
xexp[-S(xi,Xi,w,Q,tji,tj2)] , 



(57) 
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where 



H2 = -^as{t2l)as{t22)Ko y^/xlx[Qblj Kq yy/x^Qb2 

are derived from the Fourier transform of Eq. (^7|), and the expressions for "ifj are the same as those in Eq. ( ^8|) 
but with il'{ki, P, /i) replaced by (f>{xi, w). The variable 9 is the angle between bi and b2. The Sudakov exponent 
S is written as 



S{x^,x[,w,Q,tji,tj2) = ^^s{xiQ,w) + 3 

1=1 

3 

■^s(xj(5,w) +3 



ill 



!=1 



7? 



(58) 



Similarly, the two-loop running coupling constant in Eq. (19) is employed to derive the complete Sudakov 
factor. This is an improvement of the analyses in [^6[ 
The evolution of /jv is given by 



fNiw) = /at (mo) 



asiw) 



-I l/(6/3o) 



(59) 



with /io ~ 1 GeV and the initial condition /jv(/io) = (5.2 ± 0.3) x 10 ^ GeV^ For the wave function (j), we 
consider both the CZ model 12^] and KS model EH) . They are decomposed into 



ix,,w)^(b^^'{x,)J2Nj 
3=0 



asiw) 



as(Mo) 



ajAj{xi) , 



(60) 



where the constants Nj, Uj and bj, and the first six Appel polynomials Aj{xi) are listed in Table 2. The 

function (/)'^'^(a:i) = 120a;iX2a;3 is the asymptotic form of (p mentioned before. 



5.3 Numerical Results 

If we employ the Sudakov exponent in Eq. ( p8| ) directly, we obtain the results Q'^Ff — 0.55 GeV* for = 10 
GeV^ and Q'^Ff = 0.48 GeV^ for = 40 GeV^ by using the KS wave function, which are only about half 
of the experimental data. These values are close to those derived in j48|. However, we argue that there is a 
freedom to choose the scales for the Sudakov evolution [|lo). For example, the ultraviolet cutoff cS^Q/V^ for the 
integral in Eq. (|l7|), with c a constant of order unity, is as appropriate as ^Q/\/2. Hence, we may determine 
c by fitting the predictions for the pion form factor to the data Q^F^ ~ 0.35 GeV^ for ~ 5-10 GeV^, and 
apply it to the evaluation of the proton form factor. This procedure leads to c w 0.6. The change from c = 1.0 
to c = 0.6 enhances the results of the pion form factor by about 15%. We expect that this change in c will give 
a larger enhancement in the proton case, since the proton form factor is more sensitive to the evolution scales. 
We emphasize that the parameter c — 0.6 is introduced in order to explain the data. It is certainly legitimate 
to question this smaller value of c, and argue that higher Fock state contributions may be important. 

In the numerical analysis below we adopt the Sudakov exponent S{xi, x[,w, 0.6Q, iji, ^^2)- We investigate the 
sensitivity of our predictions to the large 61-62 region, choosing A — 0.2 GeV. If the perturbative contributions 
dominate, an essential amount of Q'^Ff will be quickly accumulated below a common cutoff be of 61 and 62. The 
numerical outcomes are displayed in Fig. 9. All the curves, showing the dependence of Q^Ff on 6c, increase 
from the origin and reach their full height at be = 1/A. A quantitative measure of self-consistency is that half 
of the contribution to Q'^Ff is accumulated in the region characterized by as/iT < 0.5. Based on this standard, 
the results above ~ 20 GeV^ are reliable. At a higher energy = 40 GeV^, the contribution from the 
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perturbative region amounts approximately to 65% of the entire contribution. It implies that the applicability 
of PQCD to the proton form factor at currently accessible energy scale ~ 30 GeV'^ may be justified. 

Note that the saturation of the proton form factor in the large h region is slower than that in the pion case, 
because the Sudakov suppression for the former is more moderate compared to that for the latter. The proton 
form factor involves more 7^, which enhance the nonperturbative contributions, and the partons are softer, 
thereby weakening the Sudakov suppression. The slower transition of the proton form factor to PQCD as 
increases can also be understood as follows. The hard gluon exchange diagrams give leading contributions of 
order (as/7r)^, rather than a^/Tr as in the pion case, so that soft contributions, always suppressed by a power 
of M^/Q^ with M a typical hadronic mass, become relatively stronger 

In Fig. 10 we compare our predictions for Q^Ff{Q^) to the experimental data, which are extracted from 
the data of elastic electron-proton cross section l27[ ^ . The match between the curve corresponding to the KS 
wave function and the data is obvious. Below 20 GeV^, the deviation of the curve from the data becomes 

larger. The CZ model generates results which are smaller than those from the KS model by a factor 2/3. The 
above analysis implies that the KS wave function is more phenomenologically appropriate. We shall adopt this 
model for the proton wave function in the future studies. The sensitivity of our predictions to the variation of 
the QCD scale parameter A is also examined. We find that the decrease of A to 0.1 GeV enhances Q^Ff by 
about 10%. 

6 Sum Rules for Compton Scattering 

In this section we review the derivation of the local duality sum rule for pion Compton scattering. The theoretical 
motivations and justifications of the validity of sum rules for the invariant amplitudes of RCS or VCS can be 
found in 

As emphasized by Nesterenko and Radyushkin in their work on the sum rule analysis of the pion form 
factor |l^ , the dominant contribution to the intermediate energy behaviour of the pion form factor comes 
from a "soft" diagram, which is nonfactorizable and is interpreted as an overlap of initial and final state wave 
functions. According to this viewpoint, the inclusion of a hard scattering in the intermediate state, as suggested 
by the PQCD factorization formalism, does not describe the correct dynamics. There is at the moment no 
general consensus on this issue, and the experimental input is going to be crucial in order to rule out some of 
the alternatives. 

We remind briefly that the sum rule analysis of VCS requires a detailed study of the properties of analyticity 
of the diagrams appearing in the Operator Product Expansion of a 4-current correlator, and this task has been 
performed in detail in The existence of a finite contour dispersion relation, which allows the application 
of the usual OPE technique, has also been shown to hold with sufficient accuracy. The sum rules obtained by 
this method, which are similar to the form factor sum rules, are now characterized by both an "s" and a 
dependence, and are shown to hold at a fixed value of the ratio —t/s. The basic idea of this approach is to 
connect the invariant amplitudes of a "true" {i.e., with on-shell pions or protons) scattering to an OPE series, 
which is expressed in terms of a perturbative "spectral function" and of the corresponding condensates. 

We show that the sum rule can be organized in a rather compact form in the Breit frame of the pion, 
where the spectral density becomes polynomial in a variable called Q. In the sum rule description of both form 
factors and Compton scattering, Q is a function of various kinematic invariants, such as the pion virtualities 
(here denoted as si, S2), plus the usual s and t invariants. In fact, we recall that the OPE side of the sum 
rule is calculated with off-shell pions, since it has to correspond to a short-distance factorization. This variable 
Q, which parametrizes the large light-cone components of the incoming (or outgoing) hadron in the scattering 
process, is used to analyze the power suppression of the various contributions to the spectral function. Q is 
then averaged in the dispersion integral over the virtualities of the hadron momenta. 

We also evaluate the power corrections proportional to the quark and gluon condensates, and give the 
complete expressions of the spectral functions for the helicity amplitudes involved in pion Compton scattering. 
A new Borel sum rule is derived and compared to the corresponding one for the pion form factor. The Borel 
transform is defined in its "contour form", and is used to reduce the "noise" or "background contribution" 
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to the sum rule. Ultimately, the success of any sum rule program relies on the choice of a "good" correlator 
with minimal noise and maximum projection onto the physical states that we intend to study. Assessing the 
error in this unobvious re-organization of the perturbative expansion, via the OPE, is nontrivial and requires 
a case-by-case study. We remind that sum rules are just a different model, closer to a PQCD description, and 
like all models, should be handled with caution. In some situations they can help us describe the underlying 
dynamics better than other models, and, as such, can be a very useful tool of analysis. 



6.1 Overview of the Method 

We start our discussion by quoting the main results of Consider the following 4-point amplitude as graphi- 
cally described in Fig. 11, 

'■{OlT{r,MJi,MM'>)lli^)) m • (61) 

where 

2 1 - 

Jfj, = -^ujf_,u - -d'jf^d , 77^ = ujs-fad , (62) 

are the electromagnetic and axial currents, respectively, of the u and d quarks. The choice of this specific 
correlator is motivated by the particular process we are considering, namely, pion Compton scattering. The 
pions of momenta pi and p2 are assumed to be off-shell, while the photons of momenta qi and 52 sue on-shell. 
The invariant amplitudes Hi and H2 of pion Compton scattering are extracted from the matrix element. 



Mx^{pi,P2,qi)^i J dSe-'^'-y{p2\T{My)Jxm\pi) , (63) 
by expressing it as [|j 

M^^(pi,P2,9i) =i?i(s,i)e(i)^e(i)^ + ff2(s,0e''^^e(2)'' , (64) 

where the polarization vectors e*^^^ and e*-^-* will be defined below. appears as the residue of the double 

pion pole at Pi = p| = Tn^ in the 4-point Green function T^x^^-- It is convenient to define a double discontinuity 
[D A^x,ia of Ti^amct across the pj and pj cuts by 

a,am, = r++-r+_-r_+ + r__ , (65) 

where we define r+_|_ = Ti^XficriPi + *£;P2 + *^)- 

The poles may be isolated from the Xq — > — oo, yo oo limit of T^xtj.a in Eq- (|6l|), by inserting complete 
sets of states in the resulting time ordering: 



r.A;.. = j d^xd^yd^zd^p\d^p'2e{yQ)e{~XQ) 

■^^-ipi-=c-iqi-z+ip2-y 

xS{p'l-ml)6{p'2'-ml){0My)\p'2) 

x{p'2\T{Mz)Jx{0j)\p[){p[\vi{xm + r™"*, (66) 

where F™"* includes contributions from higher states and other time orderings, that do not contribute to the 
double pion pole. We express the matrix elements of the axial currents as 

(0h.(2/)b2>=^AP2.e-^P'^■^ (67) 
(pihi(^)|0)=-^^p'l^e^^'^■^ (68) 
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which relate the single pion states to the vacuum at times ±00. Wc then derive the relation 

A^A^a = f^pif,p2a{2TTfS{pl-ml)S{pl-ml) 

X M,A(pi,P2,gi) + A™"* , (69) 

with A'^°"'* the double discontinuity of the continuum contribution F'^""*. 

On one hand, the correlator in Eq. ( |6l] ) is expanded in terms of the diagrams in Figs. 12. Specifically, Fig. 12 
represents the perturbative contribution to the expansion, and Fig. 13 describe the power corrections, which 
are proportional to the condensates of quarks and gluons, respectively. The power corrections parametrize 
the nonperturbative structure of the QCD vacuum. When the two external photons are physical and all the 
Mandelstam invariants are moderately large (see also the discussion in Q), the contribution to the leading 
scalar spectral function comes from the diagrams in Fig. 12. These diagram, compared to the corresponding 
one in the form factor case, contain an additional internal line (the line that connects the two photons), which 
is far off-shell. We can confirm that the process is, indeed, of hand-bag type when s and t are moderately 
large, with s ^ t. As discussed in this picture breaks down at small t, near the forward region. Note 
that Fig. 12(a) does not contain any exchanged gluons, differently from the expansion based on the PQCD 
factorization theorem where the one-gluon-exchange diagrams are supposed to dominate. 

On the other hand, the resonance region is modeled, as usual, by a single resonance (the double pole at the 
squared pion mass mj) and a continuum contribution starting at virtuality sg. Information on the invariant 
amplitudes of pion Compton scattering is obtained from the residue at the double pion pole of the correlator 



in Eq. (61) 



t=iP2-pif =si+S2-2Q^-^. (71) 



with 



We choose the pion Breit frame, in which the pion momenta are parametrized as 

Pi=(3w+^w, = + (70) 

The dimensionless vectors v and v are on the light cone, obeying v"^ = v'^ — and v ■ v = 1. The momentum 
transfer t is expressed in terms of Q^, which acts as a large "parameter" in the scattering process, as 

\2 _ „ , „ 0^2 flf2 
2Q 

Notice t = —2Q^ for si = S2 ~ 0. is then given by 

^^{si+s2-t + 5) = ^{s + u + S) , (72) 

S = Vi^i + S2-tr-is,S2 = ^^\''^'''' . (73) 

In the frame specified above, all the transverse momenta are carried by the two on-shell photons with physical 
polarizations. The advantage of using such a frame is that the spectral densities are rational functions of the 
variables Q^, s, si , and S2- They are also symmetric in si and S2, as expected from the time reversal invariance 
of the scattering process. This provides a drastic simplification of the analysis. 

The matrix element in Eq. (|63[) admits various possible extensions in terms of the polarization vectors, when 
the two pions go off shell. The extension must be performed in such a way that the requirements of the physical 
polarizations for the two photons, 

e^-gi e(')-g2 = 0, 

eW-e(^') = -S,,, (74) 

are still satisfied. A suitable choice of the polarization vectors will simplify the calculations. The two vectors 
e^^-* and e^^-* that we shall use differ slightly from those given in Q: 

(1)A ^ (2)A ^ 
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with 

^ gA^.pp^^^^^^ ^ pA ^ ^^pA ^ ^^pA ^ |pA ^ (7g) 
Vi = pi ■ P2 - S2 , V2 = Pi ■ P2 - Si , (77) 

a = ^(i^ipi • i? + • i?) , (78) 

R = qi+q2 , r = q2-qi . (79) 

In the Breit frame e^^^ contains only transverse components 

eW = (0,0,e^'^), eV'^-g^-O. (80) 

At the pion pole, where Pi=P2 = 0, one gets = 1^2 = ^- With this extension, the polarization vectors e^'^ 
still satisfty the requirements in Eq. ([T^), which hold for all positive si and S2, whether or not they are equal, 
and for all vi and V2 as specified in Eq. (|77|). 

The residue at the pion pole is isolated from the exact spectral density of Eq. (^ij) from a phenomenological 
ansatz of the form "1 pole plus continuum", 

Ap,.A<T = flpit,P2a {2-Kf6{p\ - ml)S{pl - ml) 

xAf,A(pi,P2,gi) + AP"* +h.r. (81) 



where M^x has been defined in Eq. (|6^), Ap°''* is the double discontinuity of the multiparticle states of the 
correlator in Eq. (|6l|), and h. r. denotes the contributions from higher resonances. In the derivation of the sum 
rule it is assumed that other single particle states, corresponding to the double poles located at momenta Pi 
and P2 equal to the masses of higher resonances, give negligible contributions. Therefore, A^j^ao- contains the 
most valuable nonperturbative information on the invariant amplitudes Hi and H2 at the pion pole, although 
mixed with the contributions of the higher states. 

A suitable way to remove the contributions of single particle states from Eq. (^l]), but still keep the pionic 
contributions {Hi and H2), is to use a projectior n, which satisfies the requirements Q|, 

n ■ qi — n ■ {pi ~ P2) = = . (82) 

A possible choice is 

3(2)±ie(i)y\ (83) 
The projector = (n+,n~,nx) then has the components 
= (viPl^ + V2P2^ + ^{qi^ +q2^)) . 

Note that only the transverse components of n'^ are imaginary. 

For pion Compton scattering, the two projected spectral densities are then given by the expressions, 

A(i2)(p2^s,i) = f^n-pin-p2{2TT)H{pl-ml)S{pl-ml) 

X [Hi{s,t)+H2is,t)] +A(i2)p-t ^ (85) 
AW(p2,s,i) = f^n-pin-p2{2Trf5{pl-ml)6{pl~ml) 

xi?i(s,i) + A^i)?"'-* . (86) 
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6.2 Spectral Functions from Local Duality 

We present the complete evaluation of the lowest-order perturbative spectral functions A'^^^p"''* and A'^^^^p"''*, 
which appear in the sum rules for Hi and Hi + H2, respectively. In a more complete analysis of pion Compton 
scattering an asymptotic expansion of the lowest-order spectral functions is no longer sufficient. While the 
dominant 0{1/Q*) and the next-to-leading contributions can be obtained by a simple power counting of the 
various integrals, the 0{1/Q^) contributions require a full-scale calculation of all the terms from Fig. 12. This 
turns out to be a rather lengthy task. 

The double cut contribution leading spectral functions A^^^p'"'* and A^^^^p"''* are calculated in terms of the 
3-cut integrals 



I[f{k',k-pi,...)] = / d^kf{k^k-pi,...) 



S+{k^)S+{{pi~kr)6+{{p2-kr) 



iPi 



91 J 



(87) 



/'[/(fc2, k ■ pi, ...)] = / d%f{k', k ■ pi, ...)S+{k')S+{{pi - kf)S+{{p2 - kf) . 



The analogy to the form factor case is clearly due to the fact that Eq. (|7|) approximates the total discontinuity 
along the positive si and S2 axes by a "3-particle cut" , with the top line of Fig. 12 kept off-shell. The dominance 
of this discontinuity compared to others can be easily proved. Other cuts are either subleading, relevant to 
forward scattering -not considered in this approach-, or not in the physical region. The evaluation of the 
integrals in Eqs. (^^ and (|88|), though not trivial, are simplified in the Breit frame. 
The final expression for A^^^^p®''* is given by 



^(12)port ^ 



3(27r) 



{AQ^ - siS2f{2Q^s - S1S2) 
, %Q''{u-2Q'')Ri2{Q'',u,si,S2) 



{AQ^ - siS2f{2Q-^u- S1S2) 



with 



Rl2iQ^, S, Si, S2) = ^ Kis, Si, S2)Q^" 



(89) 



(90) 



a polynomial of Q^. i?i2(Q^, u, si, S2) in the second term in the brackets has a similar expression with the 
argument s in Eq. ( pC| ) replaced by u. This term corresponds to the crossed diagram Fig. 12(b), where the two 
photon lines are interchanged. The explicit expressions of the coefficients 6„ can be found in ||7|. 
Using the relation 



(n ■ Pi)in ■ P2) 



s{s + t) 



(91) 



(92) 



at the double pion pole, the expression of the sum rule becomes 

f^[Hi{s,t)+H2{s,t)]= r dsi r ds2P^'^^''''-\si,S2,s,t) , 

Jo Jo 

with the spectral density 

^(12)pert ^ _L^A(l^)^'-*(si,S2,.,t) . (93) 
s(s + tj 

The derivation of Eq. ( |9l] ) is referred to Appendix A. By a simple power counting in Q^, it is easy to observe that 
the spectral function in Eq. (|9|) behaves as asymptotically, suppressed by one more power of i, compared 
to the leading perturbative predictions from the dimensional counting rules This behaviour, noticed in 
is similar to that of the pion form factor derived from QCD sum rules |Lq| . 
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The next point to address is the calculation of the power corrections Q. Since the top of the handbag 
diagram in Fig. 12 is far off-shell, the contributions from the condensates are similar to those in the form factor 
case. Therefore, our arguments about which soft contributions to be inserted in the OPE are based on the 
similarity between the form factor and Compton scattering. In the kinematic region with moderate s, t and 
u, and with fixed —t/s, all the operator insertions (condensates) which lower the virtuality of the quark line 
connecting the two photons are not taken into account. In calculating the power corrections, the general tensor 
operator G^^uGap is projected onto its corresponding scalar average (G^) = (G^^G"'^") by the relation 



(O|G^,G„0|O) = h^{Q\Gl„G-P-mg^^g,p - g^pg,^) . (94) 



with Ic the 3x3 unit color matrix and G^^ = G°^A°/2, A"^ being the Gell-Mann matrices. Equation (94) can 
also be recast into the form 

{0\G%Glp\0) = l^{Q\G%G'^P^mg^o.gup - g^^pg^^) , (95) 

using the properties of the Gell-Mann matrices. 

We work in the Fock-Schwinger gauge defined by 

zt'A^iz) = , (96) 

and employ a simplified momentum-space method to evaluate the Wilson coefficients |]5T| . For massive quarks, 
the fermion propagator in the external gluonic background is written as |5l|| 

i> + m , 1 , (7" ^I^Ga^p - m-f^G^P-fp) 



^^^^ p^-vn? ^ 2* [p^-m^'y 

- ■ (97) 



(p^ — TO^)^ 

The above expression simplifies considerably in the massless limit, where the third term corresponds to the 
insertion of two external gluon lines vanishes, and the second term corresponds to a single gluonic insertion on 
each fermion line survives. The gluonic power corrections are then given by 

A(i2)siuon_ 20a. Q«rgl"°"(g^ 2) 

277r2s(2Q2_si)4(2Q2_52)4 ' ^ ' 

with 

rSi"°"(Q2^s,si,S2) = 8Q^ - 24QS + 14Q2s2 _ + rgSi 

+2Q^ssi ~ Ss^si + irg^Si - 6Q^ss2 ~ s'^ s-i . (99) 

Equation ( p8| ) shows that the power corrections for the new sum rules are calculable systematically in full 
analogy to the form factor case. The power-law fall off of such contributions are 0(1/(5^), suppressed by l/Q^ 
compared to the perturbative spectral function. 

The quark power corrections can also be calculated . The average of the 4-quark operators is written as 

- ^^{{i^i^fwH^'^b^.bp^ ~ s'^^s'^s^^sps) . (100) 

Combining all possible insertions consistent with the handbag behaviour of Fig. 11 we derive 

(suS2,s,t) 27Q2553(-2Q2 + s^)2,3(_2Q2 + ,2)2 • ^'^'^ 

Since the close form of the Borel transformed quark power corrections can be obtained, as shown below, we do 
not present the explicit expression of r'^^'^'^^iQ^ , Si, S2) here. 
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Adding together Eqs. (|89|), ( ^8|) and (101), the OPE of the spectral density is summarized as 



-t 



s{s + t) 



-as{G 
2 



2, 20 Q«rgi"°"(Q2,s,si,S2) 



27^2 5(2g2_5^)4(2g2_ 52)4 

107ra,r9"'^'^'=(g2,si,S2) 



S{si~pI)S(s2-pI) 



(102) 



Equation (102), once integrated along the (0, A2) cut in each of the si and S2 planes, gives the sum rule 



J 7T 



H,{s,t) + H2{s,t) 



9 9 

P1P2 



dsi 



p<-'^\si,S2,S,t) 



(103) 



^(•Sl -P?)(S2 - pi)' 

where we have included all the contributions to the sum rule under a single integral. Certainly, we still have to 
apply the Borel transform to ( |103D in order to obtain the Borel sum rule [Q. Here A defines the radius of the 
contour for the dispersion relation as shown in Fig. 14, and also fixes the parameter of the Borel transform. 
After applying the Borel transform, Eq. (103) becomes Q 

W [Hiis,t) + H2{s,t)] ' 



dsi 



ds-. 



320g"^P"t(Q2,s,,s2) 



3s(s - 2Q2)(4Q4 _ s^s2)5(-2Q2s + S1S2) 



277r2 



Q«rg'— (Q2,,,,,,,2) ^-f.,+.,)/M- 

's(2Q2_5^)4(2Q2_52)4e 



1 - e 



-(A2-si)/Af= 



with 



Fo{M,s,t) 

80 {-AM^s - 8M^s^ - 2s3 _ 2M2st - s^t + 2MH^ + 2st^ + t^) 



(104) 



27 

FiiM,s,t) 
F2(M,s,<) 



8 (2A-/2 + S + t){Al\/Ps + 2s2 -St- t^) 



(105) 



27 MH"^ 
320 s 

M being the Borel mass. 

We remind that Borel transforms, in a finite region of analiticity, should be replaced by their representation 
as contour integrals Q . A study of the stability analysis of this sum rule shows that A can be chosen to coincide 
with the duality interval sq [p^ . We compare the above formula with the sum rule for the pion form factor, in 
which the integrals for the gluon and quark power corrections can be explicitly performed 



16 1/ 



127rM4 



dsi 



ds2P-KS.{,si,S2,t)e 



-isi+S2)/M^ 



13 M2 J 



(106) 
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with 



3 



27r4 



[(si + S2-t)2-4siS2)]l/2 



(107) 



Note that the power corrections in Eq. ( |l06|) tend to grow rather quickly with Q^, and invahdate the sum rule 
at beyond 4-5 GeV^, while the perturbative contribution decreases fast as These behaviours set a 

tight con straint on the validity of the sum rule. The asymptotic dependence of the perturbative term in 
Eq. (104) is similar to that in Eq. ( 106 ), but the power corrections are more stable. This fact is due to the extra 
angular dependence of Compton scattering. To have a more precise and definite answer to the above issue, a 
stability analysis of the sum rule in Eq. (104) and a determination of the Borel mass M are required. This 
study has been performed in [O. 



7 Perturbative Pion CS: a Comparison with QCD Sum Rules 

In this section we apply the resummed PQCD formalism to pion Compton scattering. We review the calculation 
of the invariant amplitudes Hi and H2 , and compare them to the sum rule results obtained in the previous 
section. A comparison of perturbation theory and QCD sum rule prediscutions is presented. 



7.1 Factorization Formulas 

The reasoning leading to the factorization of pion Compton scattering is similar to that for the pion form factor. 
We consider first the factorization formula for H = Hi + H2, which includes the Sudakov resummation of the 
higher-power effects at lower momentum transfers. Basic diagrams for pion Compton scattering in the PQCD 
approach are shown in Fig. 15, which differ from those in Fig. 12 by an extra exchanged gluon. The contribution 
to the hard scattering amplitude from each diagram in Fig. 15, obtained by contracting the two photon vertices 
with —g^'^, is listed in Table 3. All the contributions can be grouped into two terms using the permutation 
symmetry. We have 

2 „i 



H{s,t) = / dxidx2 4>{^i)4'{x2) I bdbHi{xi, s,t,b) 

i^-^ Jo Jo 



X exp[-S'(a;„ 5, Q)] , (108) 

similar to Eq. (|29|). Here the pions are assumed to be on-shell, pf = 0, and t = —2Q^ is obtained from Eq. ( [7l| ) 
by setting Si to zero. The Sudakov factor exp(— S") is the same as that associated with the pion form factor in 
Eq. (|30|). The function (j) is taken as the CZ model |25| in Eq. (||), since it gives the predictions for the pion 
form factor which are close to those from the wave function <f) in Eq. (|3^). In this way we simplify the analysis. 
The Fourier transformed hard scattering amplitudes Hi are given by 

167TCF(el + el)as{wi) ( r^\(\{\-xi)t^u\\{\-X2)t + u\ 



Hi 



(\-Xi)(\-X2) 

[(l-a;i)t + ,s][(l 



from the classes of Figs. 16(a)-16(c 




(109) 



H2 = S27rCFeuedasiw2) e{-r2)Ko[y/\r2\b) +0ir2)'-^H^''HV^b) 
1 (1 + X2 — XiX2)t^ + (1 + a;2 — Xi)ut 



Xi{l-Xi) X2{l-Xi)s'^ 

1 [1+ Xi— XiX2)t^ -f (1 -I- 2:1 — X2)st 

X2{l-X2) Xi{l-X2)v?' 



(110) 
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from the classes of Figs. 16(d) and 16(e), with 
ri = XiX2t, r2 = XiX2t + Xiu + X2S . 



(Ill) 



Kq and hI^^ in Eqs. (109) and ( |110| ) are the Bessel functions in the standard notation. The imaginary contri- 
bution appears, because the exchanged gluons in Figs. 16(d) and 16(e) may go to the mass sheh (r2 = 0). The 
argument wi of Us is defined by the largest mass scale in the hard scattering. 



Wi 




W2 



(112) 



As long as h is small, soft r; does not lead to large as. Therefore, the nonperturbative region in the modified 
factorization is characterized by large h, where Sudakov suppression is strong. Equation ( |108| ), as a perturba- 
tive expression, is thus relatively self-consistent compared to the standard factorization. Since the singularity 
associated with r2 = is not even suppressed by the pion wave function, Sudakov effects are more crucial in 
Compton scattering than in the case of form factors. 

Following the similar procedures, we derive the first invariant amplitude Hi. The extraction of Hi can be 
performed by contracting the two photon vertices with e'^^^^e'^^^'^ . The derivation of Hi is much simpler than 
that of 772, because e^^-* is orthogonal to all of the momenta Pi and qi. The contribution to the hard scattering 
amplitude associated with Hi from each diagram in Fig. 16 is also listed in Table 3. The modified perturbative 
expression for Hi is given by 



2 „i 

i?i(s,t) = ^ / dxidx2 (j){xi)cj){x2) 
1=1 -^0 

X exp[-S'(a;i,6,Q)] , 



bdbHl^\xi,s,t,b) 



with 



^(1) ^ SnCriel + e^)as(wi) 



(l-a;i)(l-a;2) 
from the classes of Figs. 16(a)-16(c), and 



-H h4-2a;-2?/ 

■s u 



(113) 



(114) 



H. 



(1) 



167rC_Fe„erfas(w2) 
1 



e(-r2)K^W\r2\h] +e(r2)yi/^'^ (V^5) 



1 



t 



t 



_Xi{l-Xi) X2{l-X2) X2{l~Xi)s Xi{l-X2)u 

from the classes of Figs. 16(d) and 16(e). 



(115) 



7.2 Numerical Results 

Based on the sum rules and the resummed perturbative formulas for Hi and H2 of the previo us se ctions, we 
compute the magnitudes of these two helicities. Sum rule predictions arc obtained from Eq. (104) with the 
substitution of so — 0.6 and — 4 GeV^ 0]. The resummed PQCD formula for Hi in Eq. (112) is evaluated 
numerically, and H2 is derived by H2 = H — Hi. Results of Hi and H2 at different photon scattering angles 
9* are shown in Fig. 16 and 17, respectively, in which \Hi\ denotes the magnitude of Hi. Note that the sum 
rule predictions for Hi are negative, and those for H2 are positive. It is observed that the sum rule results 
decrease more rapidly with momentum transfer \t\ |Q, and have a weaker angular dependence compared to 
the PQCD ones. The PQCD predictions are always larger than those from sum rules at 0* = 50° {—t/s = 0.6), 
and are always smaller at 6* = 15° {—t/s = 0.2) in the range 4 < |t| < 16 GeV^. It implies that large-angle 
Compton scattering might be dominated by perturbative dynamics. These two approaches overlap at 9* = 40° 
{—t/s = 0.5) and at |i| = 4 GeV^, showing the transition of pion Compton scattering to PQCD. The transition 
scale is higher at smaller angles. 
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With these results for Hi and H2, we compute the cross section of pion Compton scattering. The expression 
for the differential cross section of pion Compton scattering in the Breit frame is derived in the Appendix B as 



da |i/i|2 + |i/2|Vs-?/^ ' 



dcos6'* 2567rt 



(116) 



which yields the results exhibited in Fig. 18. In the angular range we are investigating, the PQCD and sum rule 
methods predict opposite dependence on 9*: the PQCD results increase, while the sum rule results decrease, 
with the photon scattering angle. The reason for this difference is that the increase of the amplitudes with 
9* is not sufficient to overcome the increase of the incident flux (see the Appendix B). At a fixed angle, the 
differential cross section, similar to Hi and H2, drops as \t\ grows. Again, the transition scale is around 4 GeV^ 
for 9* = 40°. 

The differential cross section of pion Compton scattering from a polarized photon has been analyzed based 
on the standard factorization formula in ||54|. The relevant scattering amplitudes and cross section can be 
expressed in terms of Hi and H2 evaluated here . It was found that our predictions are consistent with those 
obtained in However, in ||5^ the coupling constant as is regarded as a phenomenological parameter and 
set to 0.3, while here we consider the running of as due to the inclusion of radiative corrections, and its cutoff 
is determined by the Sudakov suppression. Furthermore, our perturbative calculation is self-consistent in the 
sense that short-distance (small-6) contributions dominate. 

In our analysis, which is restricted to lowest order in a^, the sum rules for the helicities are real, and 
therefore, the issue related to the perturbative and nonperturbative nature of the phases of Compton scattering 
is not addressed. This issue needs a further study. 



8 Conclusion 

In this work we have recalculated the pion form factor using the PQCD factorization theorem, which modifies 
the standard exclusive theory [ p^ by including the transverse degrees of freedom of a parton and the Sudakov 
resummation. The large logarithmic corrections at the end points of the parton momentum fractions are 
summed into an exponential factor, which suppresses the nonperturbative contributions from soft partons. 
With the help of the Sudakov suppression, PQCD enlarges its range of applicability down to the energy scale 
of few GeV. Compared to the previous analyses, we take into account the evolution of the wave function and 
the more complete expression of the Sudakov factor derived from the two-loop running coupling constant. The 
experimental data are then explained. 

By extending the resummed formalism to the proton Dirac form factor and choosing an appropriate ultra- 
violet cutoff for the Sudakov evolution, we have found that perturbation theory becomes self-consistent for 
beyond 20 GeV^, and our predictions from the KS wave function match the experimental data. Compared to 
the standard leading-order and leading-power descriptions [|l^, |l9|, ^ , where the transition to the perturbative 
behaviour is relatively slow, and the predictions are not reliable even at the highest available energies our 
approach provides an improvement of the theory for QCD exclusive processes. In the present calculation we 
have modified the infrared cutoff for the Sudakov evolution, and render it the same as the factorization scale for 
the proton wave function. By this modification, the soft divergences from the factorization scale close to Aqcd 
are completely removed by the Sudakov suppression. Our analysis shows that the KS wave function is a better 
choice compared to the CZ wave function, since the latter gives results which amount only to 2/3 of the data. 

We have reviewed the previous works on the local duality sum rule for pion Compton scattering and explained 
the derivation of the leading spectral function for the two invariant amplitudes. The leading nonperturbative 
corrections to the sum rules can also be calculated in a systematic way. It has been shown that the sum rules 
for pion Compton scattering, though more complicated, has the same power-law fall-off as that of the pion form 
factor in the related sum rule analysis [|l6j . 

The application of the resummed PQCD formalism to pion Compton scattering is also reviewed. We have 
given the modified perturbative expressions for the invariant helicity amplitudes. Comparing the predictions 
for pion Compton scattering from the resummed PQCD formalism with those from the sum rule analysis, we 
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observe clear overlaps between the two approaches. The numerical results show that the transition to PQCD 
appears at \t\ = 4 GeV^ and at 9* = 40° in the Breit frame. This suggests that sum rule and PQCD methods 
are complementary tools in the description of exclusive reactions, and the power-law fall-off of the relevant 
amplitudes helps us locate the transition region. 

A more convincing justification of our conclusion will come from the study of other more complicated 
processes, such as proton Compton scattering |^. We have started a complete analysis of VCS using the 
resummed perturbation theory and QCD sum rules, and the results will be published in the near future. Finally, 
we expect that our methods will find applications in the entire class of elastic hadron-photon interactions at 
intermediate energies. 
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9 Appendix A. Kinematics 

Let qi and q2 be the momenta of the incoming and outgoing photons, respectively, which are assumed to 
be on-shell {qf = =0). Let pi and p2 be the momenta of the incoming and outgoing pions. The external 
pion states are off-shell and are characterized by the invariants si = and S2 = pI- We define the Mandelstam 
invariants 

s = {pi + qif , t = (g2 - qif , , u = (p2 - qif , (A.l) 
which satisfy the relation 

s + t + u = SI + S2. (A.2) 
It is also convenient to introduce the light-cone variables for the photon momenta: 

qi = qtv + q^v + qi±_ , 

i; = -^(l,l,Oi), ^^-L(i,_i,Oi), 

gi± • n± = , (A.3) 
with the convention 

= + v^^{v'-v^). (A.4) 

Covariant expressions for q^ are found to be 

In the Breit frame of the incoming pion we easily get 

u={p2-q,f = 2Q^-s+'^. (A.7) 
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As discussed in Sect. 6, a crucial step in the derivation of the sum rules is the appropriate choice of the 
projector n'^. Its real part is parametrized as 

nn^{n+,n-,^^). (A.8) 

We easily obtain the expressions 
{4Q^ - iQ\s + S1S2) 



)2 



2Q 

:2Q^-5i)(2( 
4Q4 



(A.9) 



(2g2 - s)(2Q2 _ si)(2g2 - S2)(2g2s - S1S2) 



p. = ^1-^ tn^ ^ . (A.io) 



Note that both and a are negative. The relation 

n-pi = (l/V'^) (n+pj; +n_p+) , (A.U) 
(s-2Q2)(-2Q^s + siS2) 

2g2 

together with the condition n ■ pi = n ■ lead to 

1 (2g2-.si)(2g2-s2) 



n-pin-p2 (2g2 - s)(2g2s - S1S2) ' 
At the pion pole (si = S2 = 0, g^ = — </2), we have 
1 -t 



n ■ Pin ■ p2 s{s + t) 
which is the result quoted in Eq. (^). 



(A.12) 



(A.13) 



10 Appendix B. The Scattering Amplitudes and Cross Section 



In this appendix we derive the expressions for the amplitudes of pion Compton scattering and the cross 
section appearing in Sect 7. The differential cross section of Compton scattering is defined by 



\M\^ 
da = ^^dg , 

where A4 is the scattering amplitude, F is the incident flux, and Q is the phase space of the final states, 

dp2 dq2 



dg = (27r)V4)(pi+gi 



P2 - q2) 



(27r)3pO (2^)3g§ 

F is defined by 

F ^ Wp, ^ Vq,\2p°2q1 = Ap ■ q = 2s , 



(B.I) 



(B.2) 



(B.3) 



with Vp ~ p/p'^ the velocity of the incoming particle. It is easy to observe that F increases with the photon 
scattering angle in the Breit frame. Combining Eqs. (B.2) and (B.3), we obtain the general expression 



da — — — dcoa9 
S2ns 



(B.4) 
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with 9 the center-of-mass scattering angle. 

Substituting cos 9 = {t — u)/s into eq. (B.4), we have the expression which is invariant in both of the center- 
of-mass and Breit frames. Then Eq. (B.4) can be easily converted into the one in the Breit frame using the 
relation sin6'*/2 = —t/{s — u) with 9* as defined before. We have 



da 



dcos 9* 



1287rt 



where the scattering amplitude M is given by 



(B.5) 



(B.6) 



with All 
Eq. (pl6|) 



fj^i, defined by Eq. (pi]), and et the polarization vector of the photon in the state T. Inserting |A^p = 
(- |iJ2h)/2 for pion Compton scattering from an unpolarized photon into the above formula, we obtain 
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Tabic 1. The expressions of the integrand Ya' /s'-y' Ha' is'YafS-yYaiS-y for the diagrams in Fig. 8. ■!/'i23 is the brief 
symbol for xl'{ki. ^2- k:>,. P. //.) and V'j^.-j for V'(k'\ ■ k!,. k'^. P'. //). 



Diagram YHY/{4TT^alf^/27) 



(a) 
(b) 
(c) 
(d) 



eu(-0123'0i23 + 471237(23) 



(g) 

(i) 



(k) 



(1 - - :ri)[(l - x,){l - x[)Q^ + (kr, - k'^J^][xsx'^Q^ + (kr^ - 

e»(V--'l23-0l23 + 4ri23ri'23) 

(1 - Xi)(l - X[)[{1 - Xi)(l - X[)Q^ + (kT, - k^J^][X2X',Q^ + (kT, - k^J2] 

4Ti23T{23 

(1 - Xi)(l - 4)[X24Q' + (kr, - k^J2][^3a./ Q2 + (k^^ _ k^J2] 

e«V'l23'^i'i23 



(1 - Xi){l - X',)[X2X'^Q^ + (kT, - k^J2][3.3a./Q2 + (kT3 - k^^)^] 

(e) 

(f) 

e«V'213V'213 



(1 - X3)(l - X'^)[X2X'^Q^ + (kT, - k^J^][X3x',Q^ + (kr, - k^J2] 

ed(V'l23V'l23 + V'213'^213) 

(1 - X3)il - x':,)[x2x',Q^ + (kT, - - a;3)(l - x',)Q^ + (kr, - k!j.^y 

edV'l23V''l23 



(1 - X3)(l - x[)[xix[Q^ + (kr, - k^J2][^24Q' + (kr, - k^J2] 
(j) 

ed4ri23Ti'23 



(1 - X2)(l - a;'i)[a;ia;iQ' + (kr^ - K^?][x2x'2Q^ + (kr, - k^J^j 
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Table 2. The constants and the Appel polynomials for the proton wave function (f>{xi, w) of the CZ and KS 
models ||, || in Eq. @. 



j 


aj(CZ) 


aj(KS) 











1.00 


1.00 


1 





1 


1 


0.410 


0.310 


21/2 


20/9 


Xi - X3 


2 


-0.550 


-0.370 


7/2 


24/9 


2 — 3(a;i + X3) 


3 


0.357 


0.630 


63/10 


32/9 


2 - 7{xi + X3) + 8{xl + xl) + AX1X3 


4 


-0.0122 


0.00333 


567/2 


40/9 


.Ti ~X3^ (4/3)(a;f - xl) 


5 


0.00106 


0.0632 


81/5 


42/9 


2 — 7(a;i + x^) + lAxiXs 












+ (14/3)(x? + ^i) 



35 



Tabic 3. The expressions of the hard scattering amplitudes H and H^^^ corresponding to the diagrams in 
Fig. 16. Here we define 

Di = xiX2t - (kTi - IctJ^ , 

D2 = XlX2t + XlU + X2S — (kTi — kTj)^ . 

Diagram H/jlQiTasCF) H'^^^ /{STrasCp) 



(a) 

(b) 

(c) 
(d) 

(e) 



-el[{l - xi)t + u][{l - X2)t + u] 



{l-Xi){l-X2)s'^Di {I - x^)(l - X2)sD^ 

«2 „2 





{l-xi)Di 















Xi{l-Xi)D2 Xi{l — Xi)D2 

eued[(l + X2- xiX2)t^ + {l + X2- x{)ut] e„edi 



X2{1 - Xi)s'^D2 X2{l-Xi)sD2 




Figure 2: Radiative corrections to the basic scattering diagrams in Fig. 1(b). 
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Figure 3: (a) The derivative dV/dlnQ in the axial gauge, (b) The 0{as) functions K and G. 
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0.40 




Figure 4: Dependence of Q^F^{Q^) on be from in Eq. ( p5[ ) for — i GeV^ (dashed line) and for — 9 
GeV^ (solid line). 
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0.4 



(a) 



evol. 

AS 

CZ 



0.3 



O 



O 



0.2 



0.1 



0.0 



0.0 2.0 4.0 6.0 8.0 10.0 

q' (GeV') 



Figure 5: Dependence of Q^Ft,{Q^) on from (p^^ (dotted line), form (p'^^ (dashed line), and from (p in 
Eq. (|5|) (solid line). 
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Figure 6: Dependence of Q^Ft^{Q^) on for the use of in Eq. ( ]35| ) derived from Eq. ( |29| ) (solid line) and 
from Eq. ( |3^ ) (dashed line). 
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(0) (d) 



Fig. 7 

Figure 7: (a)-(d) Radiative corrections to the basic diagrams for the proton Dirac form factor. 
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Fig. 8 



Figure 8: Independent lowest-order hard scattering subdiagrams for the proton Dirac form factor. The bottom 
line represents the d quark. 
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1.0 



Figure 9: Dependence of Q'^Ff{Q'^) on the cutoff be from the KS wave function for = 20 GeV^ (dotted hne), 
for Q2 = 30 GeV^ (dashed Une), and for = 40 GeV^ (soHd hne). 
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Figure 10: Dependence of Q'^F[{Q^) on from the KS wave function (sohd hne) and from the CZ wave 
function (dashed hne). The experimental data are also shown. 
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Figure 12: (a) and (b) Lowest order perturbative contribution to the spectral density for pion Compton scat- 
tering. The dashed hues are on shell. 
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Figure 14: The integration contours for the scalar amplitude in the pf plane, showing the u = threshold. 
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1.0 




Figure 16: Dependence of (a) on \t\ derived from the modified PQCD (solid lines) and from QCD sum 

rules (dashed lines) for (1) -t/s = 0.6 {9* = 50") (2) -t/s = 0.5 (6** = 40°), and (3) -t/s = 0.2 {6* = 15°). 
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0.0 I ' ' ' ' ' ' ' 1 
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Figure 17: |t||i?2| on |t| derived from the modified PQCD (solid lines) and from QCD sum rules (dashed lines) 
for (1) -t/s = 0.6 (6** = 50°) (2) -t/s = 0.5 (0* = 40°), and (3) -t/s = 0.2 (6** = 15°). 
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10"^ I . < . < . < . 1 

0.0 5.0 10.0 „ 15.0 20.0 

|t| (GeV^) 



Figure 18: Dependence of da/dcosO* on \t\ derived from the modified PQCD (solid lines) and from QCD sum 
rules (dashed lines) for (1) —t/s = 0.6 (2) —t/s = 0.5, and (3) —t/s = 0.2. Note that the curve (1) from sum 
rules is shown by a long-dashed line. 



